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RADIATION  FROM  AN  ELECTROMAGNETIC  SOURCE  IN  A  HALF -SPACE  OF 
COMPRESSIBLE  PLASMA -SURFACE  WAVES 

by 

S.  R.  .Seshadri 
Gordon  McKay  Laboratory 
Harvard  University 
Cambridge,  Massachusetts 


^  ABSTRACT 

The  radiation  characteristics  of  a  line  source  of  magnetic  current  are 
studied  for  the  case  in  which  the  source  is  situated  in  a  half-space  of  isotropic, 
compressible  plasma  which  is  bounded  on  one  side  by  a  perfectly  conducting,  rigid 
plane  screen.  In  addition  to  the  electromagnetic  and  plasma  space  waves,  the  line 
source  excites  a  boundary  wave,  This  boundary  v/ave  is  a  coupled  wave.  It  has 
associated  with  it  both  a  magnetic  field  component  and  the  pressure  term.  This 
is  in  contrast  to  the  space  waves  which,  can  be  decomposed  into  an  electromagnetic 
(EM)  mode  with  no  pressure  term  and  a  plasma  (P)  mode  with  no  magnetic  field 
associated  with  it.  The  characteristics  of  this  boundary  wave  are  evaluated. 

The  boundary  wave  propagates  for  all  frequencies  and  the  power  carried  by  the 
boundary  wave  becomes  smaller  as  the  frequency  is  increased. 


TR396 


INTRODUCTION 

The  study  of  the  radiation  characteristics  of  localized  electromagnetic 
sources  in  ionized  gaseous  media,  known  generally  as  plasmas,  has  application 
to  the  problem  of  radio  communication  with  missiles  and  space  vehicles  passing 
through  the  ionized  regions  in  space.  In  this  paper,  the  radiation  characteristics 
of  a  line  source  of  magnetic  current  situated  in  a  compressible  plasma  medium 
of  semi-infinite  extent  and  bounded  on  one  side  by  a  perfectly  conducting,  rigid 
planar  screen,  is  investigated.  The  special  case  for  which  the  line  source  is  on 
the  screen  becomes  equivalent  to  the  problem  of  radiation  into  a  plasma  half¬ 
space  from  an  infinitely  long  slot  radiator  whose  width  is  very  small  in  comparison 
to  the  wavelength.  In  a  previous  paper  [  1  ]  ,  the  same  problem  was  investigated 
for  the  case  of  an  incompressible  plasma  in  which  the  longitudinal  plasma  waves 
cannot  be  sustained.  In  this  investigation,  the  previous  analysis. is  extended  to 
the  case  of  a  compressible  plasma  which  is  capable  of  supporting  both  the 
transverse  electromagnetic  waves  and  the  longitudinal  plasma  waves. 

Only  recently,  the  radiation  from  sources  in  a  compressible  plasma  has 
been  studied.  Hessel  and  Shmoys  [2]  have  treated  the  problem  of  radiation  from 
a  point  source  of  electric  current  in  a  homogeneous  isotropic  plasma.  Hessel, 
Marcuvitz  and  Shmoys  [3j  have  investigated  the  radiation  from  a  line  source  of 
magnetic  current  in  a  vacuum  in  the  presence  of  a  plasma  half-space.  In  this 
paper  the  radiation  characteristics  of  a.  line  source  of  magnetic  current  are 
investigated  for  the  case  in  which  the  source  is  located  in  a  half-space  of 
isotropic,  compressible  plasma  which  is  bounded  on  one  side  by  a  perfectly 
conducting,  rigid  plane  screen.  In  addition  to  the  EM  and  P  space  waves,  a 
boundary  wave  is  found  to  exist  along  the  screen.  This  boundary  wave  is  a 
coupled  wave;  it  has  associated  with  it  both  a  magnetic  field  and  a  pressure 
term.  This  is  in  contrast  to  the  space  waves  which  can  be  decomposed  into 
two  distinct,  uncoupled  modes,  namely  the  transverse  electromagnetic  (EM) 
mode  and  the  longitudinal  plasma  (P)  mode.  The  characteristics  of  this 
boundary  wave  are  evaluated. 
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The  space  wave  parts  of  the  electromagnetic  and  the  plasma  modes  do  not 
propagate  for  frequencies  less  than  the  plasma  frequency.  But  the  boundary 
wave  propagates  for  all  frequencies.  The  power  carried  by  the  boundary  wave  is 
decreased  as  the  frequency  is  increased, 

FORMULATION  OF  THE  PROBLEM 

Consider  a  perfectly  conducting  and  rigid,  planar  screen  of  infinite  extent. 
Let  a  right-handed  rectangular  coordinate  system  x,  y  and  z  be  chosen  such 
that  the  screen  occupies  the  plane  z  =  0  .  (Fig.  1).  The  half-space  z  >  0  is 
filled  uniformly  with  a  homogeneous,  electron  plasma.  A  line  source  of  magnetic 
current  is  located  in  the  plasma  at  x  =  0  ,  z  =  d  ;  it  is  parallel  to  the  y-axis 
and  may  be  represented  as 

J  =yJ  =yJ6(x)6(z-d)  (1) 

m  1  m  1  o  ' 

It  is  desired  to  investigate  the  radiation  characteristics  of  this  line  source  with 

particular  reference  to  the  boundary  waves  that  propagate  along  the  screen. 

Attention  is  given  only  to  the  steady-state  problem.  The  current  source 
is  assumed  to  have  the  harmonic  time  dependence  of  the  form  .  The 

frequency  of  the  source  is  assumed  to  be  sufficiently  high  so  that  the  motion  of 
the  ions  may  be  neglected.  In  addition,  the  strength  of  the  source  is  assumed 
to  be  so  weak  that  only  small  amplitude  waves  are  excited,  thus  justifying  the 
use  of  a  linearized  plasma  theory  [4]  .  According  to  the  linearized  theory, 
all  the  field  components  will  have  the  same  harmonic  time  dependence  as  that 
of  the  source,  namely  e  1Wt,  which  may,  therefore,  be  conveniently  suppressed. 
The  collisions  between  electrons  and  other  particles  are  neglected  and  the 
drift  velocity  of  the  electrons  is  assumed  to  be  zero. 

Let  Nq  be  the  average  number  density  of  electrons,  p  the  pressure 
deviation  of  the  electrons  from  the  mean  and  V  the  velocity  of  electrons.  Let 
E  and  H  be  the  alternating  electric  and  magnetic  fields, respectively.  The 
linearized,  time -harmonic  hydrodynamic  equation  of  motion  of  the  electrons  is 

-ium  N  V  =  NeE  -  Vp 
o  o 


(2) 


FIG.  1  GEOMETRY  OF  THE  PROBLEM 


TR396  -3- 

where  e  is  the  charge  and  m  is  the  mass  of  an  electron.  The  linearized 
equation  of  continuity  together  with  the  equation  of  state  gives 

V  •  V  =  iup  (3) 


2  .T 

a  m  N 


where  a  is  the  velocity  of  sound  in  the  electron  gas.  In  addition,  the  electric 
and  magnetic  fields  satisfy,  in  the  half-space  z  >  0  ,  the  following  time-harmonic 
Maxwell’s  equations: 


V  x  E 


iu)  A*  H  -  J 
o  m 


V  x  H  =  -iwe  E  +  N  e'V 
o  o 


(4) 

(5) 


where  u  and 
o 


£o  are  respectively  the  permeability  and  the  dielectric  constant 

of  free  space. 

The  source  and  the  geometry  of  the  problem  are  independent  of  the  y 
coordinate  and,  hence,  all  the  field  quantities  also  are  independent  of  y  .  On 


substituting 


_8_ 

8y 


=  0  in  (4)  and  (5)  ,  it  is  found  that  the  electromagnetic  field 


is  separable  into  E  and  H  modes.  Since  only  a  line  source  of  magnetic 
current  is  present  the  H  mode  is  not  excited,  and,  hence,  E  =  H  =  H  =  0  . 
From  (2),  it  is  seen  that  =  0  also  .  Only  a  single  component  of  the  magnetic 
field,  namely,  is  present. 


By  making  use  of  (2)  and  (5),  E^ , 


in  terms  of  H  and  p  as  follows: 

y  * 


E 


KdS  S, 
o  1 


4“  H 

Sz  y 


V  and 
x 


(i-e. 


V  are  easily  expressed 


N  e£, 
o  1 


9.x 


(6) 


E 


y 


y  = 


1 

_JLh  - 

(i-Ef) 

8P 

iwe  e, 

0  1 

Sx  y 

1 

O  I 

9z 

1 

e  1 

_  pr 

i(JE 

0 

8 

2 

W  m£  e,  J 

0  I 

8z  y 

N  e 
o 

9x 

e 

-  H  - 

8x  y 

iioe 

0 

a 

2  „  „ 

N  e 

8z 

•w  nie  E, 

0  1 

o 

(7) 


(8) 


(9) 
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where 


£1  =  1 


N  e 
o 


and  the  plasma  frequency  is  , 
(1)  in  (4)  yields 


(10) 


The  substitution  of  (6),  (7),  and 


where, 


+  k_ 


8x 


8z 


H  '(x,  z  )  =  -  iws0£1  Jq5  (x).  6  (  Z  -  d) 
2  2  w  2 

Vd 


~r 

c 


(id 

(12) 


and  c  is  the  velocity  of  electromagnetic  waves  in  free-space.  In  a  similar 
fashion,  the  use  of  (8)  and  (9)  in  (3)  gives 


r  8 


where 


8x 


8 


8z 


+  k 


2 

W  r 

T  el 
a 


2  l 


p (x,  z  )  =  0 

>  .  2 


« --S-) 


(13) 


(14) 


It  is  advantageous  to  rearrange  the  field  quantities  given  in  (6)  -  (9)  into 
two  distinct  groups.  The  magnetic  field  component  and  the  first  terms  in 
(6)  -  (9)  constitute  the  first  group  which  is  called  the  electromagnetic  (EM) 
mode.  The  pressure  p  and  the  second  terms  in  (6)  -  (9)  constitute  the  second 
group  which  is  called  the  plasma  [P]  mode.  The  EM  mode  and  the  P  mode  are 
written  down  separately  as  follows; 

EM  mode 


H  ,  E 
y  ’  e 


uw  £  £ , 
o  1 


~  H 
8z  y 


A 

z 


a 

ax 


*v] 


y  = 


2 

w  m£  £ , 
o  l 


s  *  ®  H 

-| —  H  -  z  — J 
8z  y  8x 


A 

X 


(15) 
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P  mode 


(i  -  ex) 


N  e£ 
o 


1 


ax 


i 

iumN  £ . 
o  1 


x 


ax 


+  J  S-] 
+ J  -H-] 


(16) 


The  subscripts  e  and  p  denote  the  EM  and  the  P  modes,  respectively.  For 

the  EM  mode,  it  is  seen  that  VxE  4  0  ,  V  •  E  =  0  and  V  •  V  =  0  .  Therefore, 

it.  is  obvious  that  the  EM  mode  has  a  magnetic  field  component  associated  with  it 

but  has  no  charge  accumulation.  Moreover,  the  EM  mode  has  a  y  component  of 

the  magnetic  field  and  this  is  perpendicular  to  the  xz -plane  in  which  the  propagation 

takes  place.  Therefore,  the  EM  mode  is  a  transverse  mode.  For  the  P  mode, 

it  is  found  that  VxE  =0,  V-  E  /  0  and  Vx  V  =  0  .  The  P  mode  obviously 

.P  P  P 

has  charge  accumulation  associated  with  it,  but  no  magnetic  field.  Since  it  does 
not  contain  any  y  component  of  the  field,  it  is  a  longitudinal  mode.  Using  the 
generalized  Poynting  vector  for  a  compressible  plasma  given  in  [2]  ,  it  can  be 
shown  that. in  an  unbounded  plasma,  the  EM  and  the  P  modes  are  uncoupled  in 
the  sense  that  the  total  power  can  be  obtained  as  the  sum  of  the  powers  in  the 
individual  modes.  However,  as  will  be  seen  later,  if  a  surface  wave  exists  along 
a  plasma  boundary,  it  is  a  coupled  wave  . 

It  is  clear  from  (13)  that  in  an  unbounded  plasma,  the  line  source  (1)  does 
not  excite  p(x,  z)  .  Therefore,  it  is  evident  from  (16)  that  the  P  mode  is 
absent  and  only  the  transverse  EM  mode  is  excited.  However,  in  a  plasma 
half-space,  on  account  of  the  presence  of  the  boundary  at  z  =  0  ,  the  longitudinal 
P  mode  is  excited.  Since  the  bounding  surface  z  =  0  ,  is  a  perfectly  conducting 
and  rigid  screen,  the  following  boundary  conditions  are  to  be  satisfied: 

Ex  (x,  0  )  =  0  (17) 


V  (x,  0  )  =  0 

z  ’ 


(18) 
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INTEGRAL  EXPRESSIONS  FOR  THE  FIELDS 

The  solution  of  (11)  and  (13)  subject  to  the  boundary  conditions  (17)  and 
(18)  will  yield  Hy(x,  z)  and  p  (x  ,  z)  .  Once  Hy(x,  z)  and  p(x,  z)  are  known, 
the  other  field  quantities  E^(x,  z  ),  Ez(x,  z  ),  V  (x,  z)and  V  (x,  z )  are  easily  obtained 
with  the  help  of  (6)  -  (9)  .  Integral  expressions  for  Hy(x,  z)  and  p(x,  z)  can  be 
easily  obtained  in  the  following  manner. 

The  geometry  of  the  problem  suggests  the  following  representations  for 
the  field  components: 

oo 

Hy(x’z)  =  !F  I  Hy(?,z)eiSx  d?  (19a) 

-OO 


00 

p  (x>z)  =  ir  I  p  (?>z> el?x 


(19b) 


-00 


Integral  representations  similar  to  (17)  are  assumed  also  for  Ex(x,  z)  and 
V  (x,  z)  .  Then  it  follows  from  (6),  (9),  (11)  and  (13)  that 


1 


E  (f ,  z)  =  t 

X  b  ’  iU)8  e ,  8z  y 

:  O  1  ’ 


8  “  W-e,  > 

*>-  -ITT? -  P(C,  ■) 


N  eE, 
o  1 


(20) 


w  me  e , 
o  1 


iwe 

-i?Hy(?,z)-  ° 


N  e  8z  P  z* 


(21) 


dz 


and 


+  s. 


H  (?,  z)  =  -  iue^j  Jo  5  (z  -  d) 


(22) 


where 


+  e  2 
2  + 


dz 


P  (t,  z)  =  0 


?i  =  +  A.2  -  t2 

f.  =  +  i  7  r2  -  k.2 


k.  >  P 

l  * 


S  >  k. 


for  i  =  e,  p 


(23) 


(24) 
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The  solution  of  (22) 

gives 
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f 

i?ez. 

A  e 

-x5  z 

+  Be 

Hy(?,z)  =  ^ 

i?ez 
Ce  e 

+  De  e 

and 

-4-  n  '(?,  d  +  .)  -  -4-  R  ( !T,  d 

dz  y  51  dz  y  *  * 


z  >  d 


z  <  d 


=  *  i  «C  E,  J 
o  1  o 


(25) 


(26) 


Similarly,  the  solution  of  (23)  gives 


if  z  -if  z 

p  (?,  z)  =  E  e  p  +  Fe  p 


(27) 


The  usual  phase  radiation  condition  is  found  to  be  applicable  for  both 

modes.  An  examination  of  the  dispersion  curves  u-k  and  u-k  for  the  EM 

e  p 

and  the  P  modes  shows  that  there  are  no  backward  wave  regions.  Therefore, 
the  radiation  condition  requires  that  H^(x,  z)  and  p(x,  z)  have  outward 
traveling  phase  fronts  for  z  tending  to  infinity.  Hence,  B  *  F  =  0  .  From  the 
boundary  conditions  (17)  and  (18),  it  is  clear  that  E^  (  J ,  0)  ,  =  f ,  0)  =  0  , 

It  follows,  therefore,  from  (20),  (21),  (25),  (27)  that 


|e(C-D)  iSd-Ej) 


WEE, 
O  I 


(28) 


and 


WE  |  E 

-  i  S  (C+D)  •  +  - -  =  0 


(29) 


The  requirement  that  the  tangential  component  of  the  magnetic  field  should 
be  continuous  at  z  =  d  gives 


ifd  i?ed  "i?ed 

Ae  e  =  Ce  +  De  * 


(30) 
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The  use  of  the  jump  conditions  (26)  in  (25)  leads  to 


A  e 


if  d 


-  (C, 


if  d 


-  De 


-if  d 


)  =  - 


u£  E, 
o  i  o 


(31) 


The  expressions  for  A,  C,  D  and  E  may  be  obtained  from  the  solution  of  the 
simultaneous  equations  (28)  -  (31)  ,  The  results  are 


A  = 


ueoeiJo 

2f. 


-i?ed  Mp-U-ej)?' 

5  + 


i£  d 
be 


Mp+(1-ei)r‘ 


B  =  0 


C  = 


■£.elJo[{e{p-(l'Cl)  t2] 

7  ]  ° 


and 


D  = 


E  = 


u  E  £,  j  i|  d 
o  1  o  *e 
-  e 


2?, 


N  e  £,  J  i£  e 
o  1  o  s 


i|  d 


I«e«p  +  <1-®1)  ?  ] 


F  =  0 


(32) 


The  substitution  of  (32)  in  (25)  and  (27)  and  the  use  of  (19)  yields  the  following 
integral  expressions  for  ( x , z )  and  p(x,  z)  : 


oo 


tt  /  \  i  r  ueoei  Jo 

Hy(x’z)  21T  J  ‘  2? 


-00 


lel  +0-5,)^  '  - 


i?x+if  z 
e  e  d? 


for  d  <  z  <  oo 


(33a) 
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Hy(x,  z)  = 


and 


00 

j_  r. 

w£  £,  J 
o  1  o 

f  -^eZ  ?eV 

(l-Cj)  S2  i«cd 

2tt  J 

-CO 

2? 

*e 

for 

0  <  z  <  d  . 

p  (x,  z)  = 

oo 

_i_  r  . 

N  e£.  J  i£ 
o  1  o 

i§ed  +  i?x  +  igpz 

2  V  J  r 

-00  [ 

I  ?  +  (1-E,)?2] 
p  e  1 

"  © 

i?x+i?  z 

e  6  d? 


d? 


(33b) 


(34) 


SINGULARITIES  OF  THE  INTEGRAND 
The  contour  for  the  integrals  in  (33)  and  (34)  is  along  the  real  axis  in 
the  f -plane  as  shown  in  Fig.  2.  It  is  desired  to  find  the  singularities  of  the 
integrands  in  (33)  and  (34).  Branch  points  are  seen  to  occur  at  ?  =  +  kg  and 
f  =  +  kp  .  With  the  help  of  (12)  and  (14),  it  is  obvious  that  the  branch  points 
at  +  kg  and  _+  kp  occur  on  the  real  or  the  imaginary  axis  of  the  f -plane 
depending  on  whether  u  is  greater  or  less  than  <jp  . 

The  poles  of  the  integrand  are  obviously  determined  by  the  roots  of 
the  equation 


(k/  -  ?2. 


(k 


s2 ) 


-r  1  /2 


+  (i  -e L)  r  =  o 


(35) 


On  eliminating  the  radical  in  (35),  £  is  found  to  satisfy  a  biquadratic  equation, 
the  solution  of  which  gives  the  following  pair  of  possible  solutions  of  (35)  as 

-  2 


'1 


=  U  +  \/U 


w 


and 


=  U 


-  \ft? 


-  w 


(36a) 

(36b) 


where 


U  = 


k  2  +  k  2 
P  e 
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2(1 - ) 
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W  = 


k  2  k  2 
— E _ 


(37) 
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ln  obtaining  (37).,  use  has  been  made  of  (10).  Only  that  root  from  (36a,  b) 
which  also  satisfies  (35)  on  the  proper  Riemann  surface,  is  the  proper  root 
of  (35).  For  finding  out  which  one  of  (36a,  b)  is  a  proper  root,  it  is 
advantageous  to  simplify  (36a,  b)  by  making  use  of  the  fact  that  the  ratio  of  the 

plasma  to  the  electromagnetic  wave  velocities  —  is  of  the  order  of  10~ 

2  c 
Hence,  — is  very  small  in  comparison  with  unity.  It  is  seen  from  (37)  that 
2  c 

U  is  very  large  compared  to  W  and,  hence,  the  following  approximation 
to  (36)  may  he  used, 


=  U  +  U  + 


Case  1 :  For  u  >  w  ,  both  k  2  and  k  2  and,  hence,  U  and  W  are 
™  "  P  ®  P 

positive  and,  hence,  the  two  roots  given  in  (38)  may  be  simplified  to  yield 


=  2U  = 


5  2  =  k  2 
2  o 


2  ,  \ 

The  root  £  given  in  (39a)  is  greater  than  k  ,  since  u>u  .  For  this 


root,  in  view  of  (24)  and  the  fact  that  k 


,  the  first  term  in  (35)  is 


negative  and  the  second  term  is  positive  and  they  can  add  to  zero.  Hence,  the 
root  (39a),  which  is  an  approximation  to  the  root  given  by  the  upper  sign  in 
(36)  is  on  the  proper  Riemann  surface.  Moreover,  this  root  is  real  and  occurs 
at  +  kg  in  the  ? -plane.  With  the  help  of  (39b),  the  root  of  ^2  is  seen  to 
be  slightly  less  than  k^  .  It  follows,  therefore,  that  on  the  proper  Riemann 
surface  defined  by  (24),  the  first  term  in  (35)  is  positive  and  so  is  the  second 


term.  They  cannot  add  to  zero  and,  hence,  the  root  (39b)  is  not  on  the  proper 
Riemann  surface.  It  is  to  be  noted  that  for  u)  >  w  ,  +  ^  is  a  real  pole  on 

the  improper  Riemann  surface. 
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2  2 
Case  2:  The  root  if,  given  in  (36a)  and  approximated  by  k  given  in  (39a) 

-  1  2  s  2 

is  real  and  positive  even  for  to  <  to  .  For  to  <  w  ,k  and  k  are  negative  and 
2  2  P  P  e  p 

if  =  k  is  positive,  the  first  term  in  (35)  is  a  negative  real  number  and  the 
s 

second  term  is  a  positive  real  number  and  they  can  add  to  zero.  Hence,  the  root 
if  j  is  a  proper  real  root  of  (35)  and  gives  a  pair  of  real  poles  in  the  lf~plane.  For 
to  <  to  ,  with  the  help  of  (12)  and  (14),  it  is  evident  that  U  >  0  and  ‘W  <  0.  Since 
W  <  0 


,  y/u“  -  W  >  U  . 


2  . 


Hence,  the  root  is  negative.  Moreover, 


^  ^ 


7  k 
2  e 


\ 

to 


<  0  and  k  ^  -  if-^  =  k  ^ 
P  -  2  p 


1 


<  0 


Hence,  the  first  term  in  (35)  is  negative  real  and  the  second  term  is  also 
negative  real  and,  therefore,  both  the  terms  in  (35)  cannot  add  to  zero.  Hence, 
for  (j  <  to  ,  +  i|  if.,  |  gives  a  pair  of  purely'  imaginary  poles  on  the  improper 

p  ^ 

Riemann  surface. 


Therefore,  it  is  seen  that  the  root  (39a)  gives  rise  to  a  pair  of  real 

poles  at  if  =  +  k  for  all  values  of  to  .  The  root  (36b)  is  on  the- improper 
~  —  s 

Riemann  surface-;  it  is  real  or  imaginary  depending  on  whether  to  is  greater 
or  less  than  u  . 


,  and  if  =  +  k  are  on 
-  P 


SURFACE  WAVE  CONTRIBUTION 

The  contour  for  the  integrals  in  (33)  and  (34)  is  along  the  real  axis  in 

the  if-plane  (Fig.  2).  Since  the  branchpoints  at  if  =  +  k^ 

the  real  axis  for  to  >  w  ,  the  contour  of  integration  has  to  be  indented 

P 

suitably  at  these  points.  As  was  mentioned  earlier,  the  phase -type  radiation 

condition  is  applicable  to  the  space-wave  parts  of  the  total  field.  Since  the 

dispersion  curves,  for  the  space  waves  for  both  the  EM  and  P  modes  do  not 

contain  any  backward  wave  regions,  the  radiation  condition  is  satisfied  by 

requiring  the  phase  fronts  of  the  EM  and  the  P  space  waves  to  be  traveling 

outward.  Hence,  the  contour  is  indented  from  above  at  If  =  -  k  and 

P 

k_  and  from  below  at  jf  =  +ko  and  +  k^  .  For  to  <  to^  ,  no  ambiguity  in 


e  •”  e  p 

the  integration  contours  at  the  branch  points  occur,  since  for  to  <  u 

kp  are  purely  imaginary. 


k  and 
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In  addition  to  the  branch  points,  the  integrands  in  (33)  and  (34)  have 

poles  at  £  =  +  k  on  the  real  axis  for  all  u  .  For  u  >  u  ,  this  pole  occurs 
s  p 

beyond  both  branchpoints.  It  is  obvious  that  the  contour  of  integration  should  be 
indented  from  above  at  £  =  -  k  and  from  below  at  f  =  k  (Fig,  2a,  b)  or 
vice  versa.  The  former  type  of  indentation  will  be  chosen  for  all  w  .  This 
ensures  the  surface  wave  to  have  outward  traveling  phase  fronts.  As  was  shown 
in  an  earlier  investigation  [5]  ,  in  a  compressible  plasma,  the  requirement  of 
outward  traveling  phase  fronts  does  not  necessarily  imply  that  the  net  power 
is  traveling  outward  at  infinityih  spite  of  the  absence  of  backward  regions  in 
the  dispersion  curve.  It  remains,  therefore,  to  be  verified  that  the  total  power 
is  indeed  traveling  outward  at  infinity. 

For  x  >  0  ,  the  integrals  (33)  and  (34)  may  be  evaluated  by  closing  the 
contour  in  the  upper  half  of  the  ?-plane.  The  contribution  to  the  integrals 
(33)  and  (34)  is  the  sum  of  the  residue  at  the  pole  ?  =  kg  and  the  two  branch-cut 
integrals.  The  values  of  the  branch-cut  integrals  depend  on  some  inverse 
power  of  x  and,  hence,  for  large  x  are  negligible  compared  to  the  contribution 
due  to  the  pole.  Thus,  for  positive  large  x  (33)  yields 


Similarly, the  contribution  from  the  pole  of  (34)  gives 


(42) 


p  (x,  z)  =  p  (x,  z)  =  P 


4^ 

s  p 
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-d/k  2-k  2 
/  s  e 


2w 


-1 

(43) 


It  is  evident  that  H  (x,  z)  and  p  (x,  z)  given  respectively  in  (40)  and  (42) 

represent  surface  waves  which  propaga.te  in  the  positive  x  direction  and  are 

exponentially  attenuated  in  the  z  direction.  The  phase  velocity  g  of  both 
s  s 

Hy  (x,  z)  and  p  (x,  z)  is  given  by 


U) 

k 


.  «  2 

1+  ~Z “ 
id 


(44) 


The  phase  velocity  is  infinite  for  u  =  0  ,  reduces  monotonically  as  the  frequency 

is  increased  and,  for  cj  =  reaches  a  value  which  is  equal  to  \flE~  times  the 

acoustic  velocity  in  the  electron  gas.  As  the  frequency  is  still  further  increased, 

the  phase  velocity  continues  to  decrease  and  asymptotically  reaches  the  value 

a  ,  the  acoustic  velocity  in  the  electron  gas,  in  the  limit  of  infinite  frequency. 

A  plot  of  phase  velocity  v  versus  the  frequency  id  is  shown  in  Fig.  3. 

s 

The  attenuation  factor  for  the  decay  of  H^(x,  z)  in  the  z  direction 

2  -1/2 

-) 


is  given  by 


k  2-k  2 
s  e 


id 

a 


1  + 


(45) 


Similarly,  the  attenuation  factor  a  for  p(x,  z)  is  obtained  as 

P  ,  -1/2 


Q  =\X 


k  2-k 2 
s  p 


i1  +  —2  ) 

'  w 


(46) 


It  follows,  therefore,  that  (x,  z)  and  p(x,  z)  which  propagate  with  the  same 
phase  velocity  in  the  x  direction  are  attenuated  differently  in  the  z -direction, 
as  seen  from  (45)  and  (46).  The  attenuation  factor  is  a  monotonically 
increasing  function  of  frequency.  It  starts  with  the  value  of  zero  for  co  =  0  , 

increases  to  the  value  —2—  for  w  =  <j  and  thereafter  continuously  increases 

a  \fl  P 
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as  w  is  increased.  However,  the  attenuation  factor  for  the  pressure  wave 
is  amonotonically  decreasing  function  of  frequency.  It  starts  with  the  value 


— E—  for  to  =  0,  decreases  to  the  value  — 2 —  for  to  =  lo  and  thereafter 
a  aN/2  P 

asymptotically  approaches  zero  as  u  is  increased  indefinitely.  A  plot  of  the 
attenuation  constants 


and  a  as  a  function  of  W  is  shown  in  Fig.  3.  For 
P 

(j  =  w  ,  both  the  attenuation  constants  are  equal.  For  u<u  .a  >  a  '  and, 
p  f  p  p  n 

therefore,  the  pressure  field  is  more  tightly  bound  to  the  guiding  surface  than 


the' magnetic  field.  But  for  u  >  , 


0.^  <  resulting  in  the  pressure  field 
being  more  loosely  bound  than  the  magnetic  field. 

It  is  desired  to  evaluate  the  total  power  radiated  per  unit  length  of  the 
source  in  the  form  of  surface  waves.  It  can  be  shown  (see  Appendix)  that  the 
total  power  radiated  by  the  line  source  is  obtained  as  the  sum  of  the  powers 
in  the  space  wave  and  the  boundary  wave  separately.  Hence,  it  is  meaningful 
to  speak  of  the  power  radiated  in  the  surface  waves  alone.  It  has  been  proved 
[  3  ,  5]  ,  that  the  time -averaged,  outward  normal  flow  of  power  through  unit 
area  is  given  by 


S  =  -i-Re[ExH*  +  pV*  ] 

The  total  power  transmitted  by  the  surface  waves  is,  therefore,  given  by 


(47) 


W  =  W, 


+  W 


(48) 


s  s 

where  W ^  and  W  are  the  net  flow  of  power  respectively  in  the  positive 

s  s 

and  the  negative  x  directions.  With  the  help  of  (47)  ,  W+  and  W  are 


easily  obtained  as 


W, 


oo 


1-  RelE’xH8’’  tpSy**  ]  dz 


00 

Re  §  [  -  EZS(X>  z>  HyS  (x,  z)  +pS(x,  z)  Vxs  ‘  (x,  z)]dz 


(49a) 
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W 


oo 

=  '  I  Re  ^  f  'EZS^X’  HyS*^X'  Z^  +  PS^X’  Z)  VXS  ^X’  dZ  > 


(49b) 


g  g 

where  H  (x,  x)  and  p  (x,  z),  for  x  >  0  ,  are  given  in  (40)  and  (42)  respectively, 


s  s 

For  x  <  0  ,  H  (x,  z)  and  p  (x,  z)  can  be  obtained  easily  as 

-ik_x  -  z  \yk  2  -  k  2 


and 


H  S  (x  z)  =  H  e  S 
y  ’  s 


P‘  (x,  z)  =  -P  e 


-ik  x  -  z 
s  ■ 


'k  2  -  k  2 
s  e 


(50a) 


(50b) 


For  x  >  0  ,  the  expressions  for  EzS(x,  z)  and  V^s  (x,  z)  are  found  by  substituting 
(40)  and  (42)  in  (7)  and  (8)  .  The  result  is 
k 


(x’z)  InTIT  HyS(x> z)  +  4-5e/P-  pS{x'z) 

o  1  o  1 


(51) 


V  s  (x,  z)  =  - 


ea, 


w  me  e, 

o  1 


H  (x,  z)  +  — 1CT - 

y  ue,N  m 

1  o 


s  s  .  . 

P  (x,  z) 


When  (40),  (42),  (51)  and  (52)  are  used  in  (49a)  ,  it  is  found  that 


W, 


H  2  + 


2 


(52) 


4ueoe1ah  ~8  T  4we1Nomap  ^s  "  2NQe  £j  HsPs. 

In  a  similar  manner  W_S  can  be  calculated  and  is  found  to  be  equal  to' W+S, 
as  expected.  Hence,  the  total  power  transmitted  by  the  surface  waves  per 
unit  width  of  the  screen  is  obtained  from  (48)  and  (53)  as 


(53) 


WS  =  2W  S  -  ^ - - - 

+  2ws  e ,  a, 
o  1  n 


s 


H  2  +  ^ - =-= - 

s  2ue  .  N  ma 


2  d-Si) 


1  o 


NoMl 


Hsps  <54» 


The  expression  (54)  for  W  can  be  considerably  simplified  when 
advantage  is  taken  of  the  fact  that  is  of  the  order  of  10-8  and,  therefore, 

c 

negligible  in  comparison  with  unity.  In  view  of  this  fact,  it  follows  from  (12), 
(14)  and  (39a)  that 


k  2  -  k  2  =  k  2  ,  k  2  -  k  2  =  k  2 
s  e  s  5  s  p  s 


u 


5 


(55) 
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The  substitution  of  (55)  in  (41)  immediately  gives 

„  _  -1 


u 


H  =  -  a  e  s ,  J 
s  o  1  o 


u 


l 


-f- 


-k  d 
s 


(56) 


In  a  similar  manner,  the  use  of  (55)  in  (43)  yields 


to 


P  =  -  N  e  J  £ , 
s  o  o  1 


u 


r~  u  4 

‘  --V 

w 


-1 


-k  d 


(57) 


With  the  help  of  (45),  (46),  (55)  -  (57)  in  (54),  it  is  possible  to  show,  after 
considerable  simplification,  that 

^ud  S2 

'  *  -&7~i 


W 


2VT 


ue  J  2  ft2  (S22  +  1) 

o  o 


(58) 


where  12  =  - 

to 

P 

• — ’  s 

It  is  seen  from  (58)  that  W  is  always  positive,  showing  that  the 
requirement  of  the  boundary  wave  to  have  outward  traveling  phase  fronts  does 
lead  to  the  fulfillment  of  the  radiation  condition,  which  requires  a  net  outflow 
of  power  at  large  distances  from  the  source. 

It  is  seen  from  (58)  that  the  power  delivered  to  the  surface  wave  is  a 
maximum  when  the  line  source  is  on  the  screen,  i.  e.  ,  when  d  =  0  .  Also, an 
increase  in  the  frequency  12  results  in  a  decrease  in  the  power  in  the  surface 
wave.  Obviously, this  analysis  is  not  valid  for  very  small  £2  ,  since  in  that 
case  the  motion  of  the  ions  cannot  be  neglected.  For  fi  <  1  ,  the  space 

waves  are  exponentially  damped  and  all  the  power  is  carried  by  the  surface  wave. 

It  is  well  known  that  in  an  unbounded  isotropic  plasma  the  EM  and 
the  P  modes  are  uncoupled  [2]  .  The  presence  of  the  boundary  provides  a 
“mechanism  of  coupling  between  these  modes.  Since  the  boundary  surface 
alone  provides  the  mechanism  of  coupling  of  power  between  the  EM  and  the 
P  modes,  it  is  intuitively  obvious  that  if  a  surface  wave  is  present,  it  will  be 
a  coupled  mode  in  the  sense  that  the  total  power  in  the  surface  wave  is  not 
equal  to  the  sum  of  the  powers  in  the  EM  and  P  modes  separately.  That  this 
actually  turns  out  be  the  case  is  indicated  by  the  presence  of  the  cross  term 
in  (54). 
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APPENDIX 

Let  E,  H,  V  and  P  be  the  total  field  quantities.  The  total  field  is  a 
solution  of  (2),  (3),  (4)  and  (5)  .  Also  let  ES  ,  HS  ,  V  and  P  be  the  field 
quantities  associated  with  the  boundary  wave  part  of  the  total  field.  The 
boundary  wave  is  a  solution  of  source-free  field  equations,  Let 

SL  =  ExHs*  +  PVS*  +  ES*  xH  +  PS*Y  (A-l) 

Taking  the  divergence  of  and  making  use  of  (2)  -  (5)  ,  it  can  easily 
be  shown  that 


-i. 

S1  = 


H 


s* 


m 


(A -2) 


The  integration  of  (A-2)  throughout  the  volume  V  enclosed  by  the  transverse 
planes  P^  and  P^  ,  the  surface  P^  at  infinity  and  the  surface  of  the  screen 
(Fig.  1)  yields 


^  I  -  Ez(x,  z)  HyS  '  (x,  z)  +  P  (x,  z)  Yj5*  (x,  z) 

P*  +P  L 
12 

-  EZS>'  (x,  z)  Hy  (x,  z)  +  PS*  (x,  z)  Vx(x,  z)J  dz 

-  I  V*  <x’z)  Jm  dv  (A-3) 

V 

The  integral  over  the  surface  Pffi  ,  which  is  located  at  infinity,  vanishes  since 
the  total  field  is  bounded  and  the  field  of  the  boundary  wave  decreases 
exponentially  in  the  z  direction.  Also, the  integral  over  the  surface  of  the 
screen  vanishes  since  z  •  S1  =  0  ,  on  the  screen,  in  view  of  the  boundary 
conditions  (17)  and  (18)  . 

The  total  field  quantities  may  be  written  in  the  form 
Hy(x,  z)  =  HyR  (x,  z)  +  HyS  (x+,  z)  ,  P  (x,  z)  =  PR(x,  z)  +  PS  (x+  ,  z) 

E  (x,  z)  =  ER(x,  z)  +  ES(x+,z),  V(x,z)  =  VR(x,  z)  +  ~VS  (x+,  z) 

for  x  >  0  (A-4) 
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and 

Hy(x,z)  =  HyR(x,z)  +  HyS(x",z),  P(x,  z)  =  PR(x,  z)  +  P®  (x“  ,  z) 

jk  r>  c  _ i  \  p  — i  s  — 

E  (x,  z)  =  E(x,  z)  +  E  (x",z),  V(x,  z)  =  V  (x,  z)  +  V  (x  ,  z) 


for  x  <  0 


(A-5) 


p  R  -»i)  ~^R 

where  Hy  (x,  z),  P  (x,  z),  E  (x,  z)  and  V  (x,  z)  denote  the  space  wave  parts 
of  the  total  field  quantities.  The  field  components  of  the  boundary  wave 
traveling  in  the  positive  and  the  negative  x  directions  are  denoted  respectively 
by  x+  and  x"  . 

In  view  of  (40),  (42)  and  (6)  -  (8)  ,  the  field  components  of  the  boundary 

wave  traveling  in  the  positive  x  direction  may  be  represented  as 

,  ik  x  ik  x 

HyS  (x  ,  z)  =  Hs(z)  e  S  ,  Ps(x  ,  z)  =  Pg(z)  e  S 

,  ik  x  ik  x 

E  (x+  .  z)  =  E  (z)  e  8  V  s  (x  ,  z)  =  V  (z)  e  S  (A-6) 

In  a  similar  way,  with  (50a,  b)  and  (6)  -  (8)  ,  it  is  evident  that  the  field 
components  of  the  boundary  wave  traveling  in  the  negative  x  direction,  are 
given  by 


H  S(x",z)  =  H  (z)  e  ^  PS(x’,z)  = 
y  s 

EzS(x',z)  =  -  Ezs(z)  e1  gX  VxS(x",z)  = 

The  substitution  of  (A-5)  -  (A-7)  in  (A-3)  , 
yields  the  following 


-ik  x 

-  Ps  (z)  e  8 

Vxs(z)e  8  (A-7) 

after  some  simplification, 


;EzR(xi,z)  Hs(z)  +  PR(x1>Z)  Vxg(z) 
-Ezs(z)  HyR(x1>Z)  +  Ps(z)  VxR(x1,zr 
-EzR(-x2-z)  Hs(z)  +  pR(-x2,z)  Vxs(z) 


ik  x. 
s  1 


o 
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+  E  (z)  H  K(-x„z)  -  P  '  (z>  VR(-x?,z) 


zs 


2’ 


ikx. 


+  2 


oo 

R«r 


dz 


-Ez»WHSW  +  P»  M  VxeM 


'-I 


H  S  (x,  z)  J  dv 
y  ’  m 


(A-8) 


v 


If  the  position  of  the  transverse  plane  is  changed,  the  second  and 
the  third  terms  on  the  left-hand  side  as  well  as  the  right-hand  side  of  (A-7)  are 
unchanged.  Hence,  it  follows  that 


ik  x 
s 


Jd2  F 


EzR(x,z)  Hg(z)  +  PR(x,z)  Vxs(z) 


-Ezs  (z)  H  R(x,z)  +  Pg  (z)  VxR(x,z) 


zs  y 

In  a  similar  manner,  it  may  be  argued  that 


=  constant  (A-9) 


-ik  x 
s 


oo 

I 


dz 


-E  R(x,  z)  H  (z)  +  PR(x,z)  V  (z) 
z  9  s  v  '  xs 


+  Ezs(z)  HyR(x,z)  +  Ps(z)  VxR(x,z) 


=  constant  (A- 10) 


The  space  wave  has  a  continuous  eigenvalue  spectrum  and  so  its  dependence  on 

ik  .x  .  t  -ik0x 

x  cannot  annul  the  factor  e  s  in  (A-9)  and  e  s  in  (A-10)  .  Furthermore, 

if  the  positions  of  the  transverse  planes  P]  and  P_  approach  respectively 

^  R  R 

x^  =  oo  and  x.,  =  oo  ,  it  is  evident  that  E^  (x,  z)  ,  P  (x,  z)j,Hy  (x,  z)  and 
YxR(x,  z)  approach  zero.  Also,  since  kg  is  real,  it  follows  that  the  constants 
in  (A-9)  and  (A-10)  are  zero.  Therefore,  from  (A-9),  (A-10),  (A-3)  and 
(A-4)  ,  it  is  evident  that 


TR396 


-21- 


=  0  (A  - 1 1 ) 

In  view  of  (A- 11)  ,  it  is  obvious  that  the  boundary  and  the  space  wave  part?  of 
the  total  field  are  orthogonal  in  the  sense  that  the  total  power  flow  is  equal  to 
the  sum  of  the  powers  separately  in  the  boundary  and  the  space  waves. 
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